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We investigate non-equilibrium two-electron transfer in a model redox system represented by a 
two-site extended Hubbard model and embedded in a dissipative environment. The influence of 
the electron-electron interactions and the coupling to a dissipative bosonic bath on the electron 
transfer is studied in different temperature regimes. At high temperatures Marcus transfer rates 
are evaluated and at low temperatures, we calculate equilibrium and non-equilibrium population 
probabilities of the donor and acceptor with the non-perturbative Numerical Renormalization Group 
approach. We obtain the non-equilibrium dynamics of the system prepared in an initial state of 
two electrons at the donor site and identify conditions under which the electron transfer involves 
one concerted two-electron step or two sequential single-electron steps. The rates of the sequential 
transfer depend non-monotonically on the difference between the inter-site and on-site Coulomb 
interaction which become renormalized in the presence of the bosonic bath. If this difference is 
much larger than the hopping matrix element, the temperature as well as the reorganization energy, 
simultaneous transfer of both electrons between donor and acceptor can be observed. 



PACS numbers: 71.27.+a, 34.70.+e, 82.39.Jn 

I. INTRODUCTION 

Electron transfer is a key process in chemistry, physics 
and biologji^i 3 -^ encountered in, e.g., chemical redox 
processes, charge transfer in semiconductors and the pri- 
mary steps of photosynthesis. In condensed polar envi- 
ronments the process involves strong coupling to the un- 
derlying nuclear motion and is usually dominated by the 
nuclear reorganization that accompanies the charge rear- 
rangement. A quantum mechanical description of elec- 
tron transfer in such a dissipative environment is given 
by the spin-boson mode&& and its variants; this model 
accounts for the essential energetics and dynamics of the 
process, such as the non-monotonic dependence of the 
transfer rate on the energy asymmetry, the energy differ- 
ence between the initial and final electronic states. 

Although standard descriptions of such processes fo- 
cus on single-electron transfer— two-electron trans- 
fer has been suggested as the dominant mechanism in 
some bioenergetic processes that occur in proteins 7 -^, 
transfer in transition metal complexe s 9 ' 10 , electrode 
reactions^, artificial photosynthesis and photoinduced 
energy- and electron-transfer processes^, biological elec- 
tron transfer chains 13 , transfer in fuel cells 1 ^ as well as 
in DNAi£. Further examples are selfexchange reactions 
such as T1(I)/T1(III) and Pt(II)/Pt(IV)^ and electron- 
pair tunneling 1 I i 1 8 1 1 9 in molecular electronic devices. 

The theoretical description of two-electron trans- 
fer dynamics differs fundamentally from its single- 
electron counterpart. More than two states have to be 
considere d 20 ' 21 and electron correlations induced by the 
Coulomb repulsion and the coupling to the environment 
need to be accounted for. Usually, the on-site Coulomb 



interaction in molecules is much larger than the inter- 
site interaction. 22 ' 23 ' 24 However, due to the polarization 
of the local environment the short range interaction may 
be strongly screened. Then, the inter-site interaction 
V can be of the same order or even exceed the on-site 
Coulomb interaction E/ , 24 ' 25 While U favors a homoge- 
neous charge distribution, the inter-site interaction V 
inclines spatially inhomogcncous charge accumulation. 
Since the non-equilibrium dynamics is governed by the 
energy difference U — V , the competition between both 
interactions influences strongly the type of charge trans- 
fer dynamics. Depending on the sign of the energy differ- 
ence a single concerted two-electron step or two sequen- 
tial single-electron steps may occur. 

In this paper we consider a system comprised of a 
donor (D) and an acceptor (A) site. They share two 
electrons which are coupled to a non-interacting bosonic 
bath. Such a donor-acceptor system has four different 
states: two doubly occupied donor (D 2 ~ A) and acceptor 
(DA 2 ~) states, and two degenerate states D~A~ with 
one electron each on the donor and acceptor site (with 
different spin). Their energy difference depends on the 
difference between on-site and inter-site Coulomb repul- 
sion as well as the bias e which we do not consider here. 
The transition D 2 ~ A — > DA 2 ~ occurs as a concerted 
transfer of two electrons or an uncorrelated sequence of 
one-electron transfer events during which the intermedi- 
ate D~ A~ is formed. The transfer rate of each electron 
may be different and shows a non-monotonic behavior 
on the energy asymmetry between the states. In this pa- 
per, we are mapping conditions under which the system 
performs concerted two-electron transfer or a sequential 
single-electron process. To this end we study the non- 
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equilibrium dynamics of the donor-acceptor system ini- 
tially prepared with two electrons at the donor site. We 
evaluate the rates for single-electron transitions and an 
electron-pair transfer in different regimes of the Coulomb 
repulsion and environmental response. 

The occurrence of such a correlated electron-pair trans- 
fer can be already understood within a donor-acceptor 
system, decoupled from the environment in which the 
strong on-site Coulomb repulsion exceeding considerably 
the inter-site repulsion. We start from a doubly occu- 
pied excited donor state D 2 ~ A compared to the D~A~ 
ground states. Energy conservation implies a concerted 
electron transfer. If the transfer matrix element A is 
much smaller than this energy difference the transfer oc- 
curs as a tunneling process of an electron pair in which 
the intermediate states D~ A~ are occupied only virtu- 
ally analogous to a "superexchange" process (see, e.g., 
Ref.H). 

In the present paper we investigate the effect of cou- 
pling to a dissipative bosonic environment, with a to- 
tal number of two electrons occupying donor and ac- 
ceptor sites. These two electrons experience the on- 
site Coulomb repulsion U when occupying the same site 
and the Coulomb repulsion V when occupying different 
sites. In this paper we restrict ourselves to the sim- 
plest case where donor and acceptor are each modeled 
by a single molecular orbital. In such a system the dif- 
ference U = U — V is crucial for the dynamics. The 
coupling to the bosonic bath has two major effects: (i) 
the renormalization 26,27 of the on-site Coulomb repul- 
sion U to U c s and (ii) dephasing as well as dissipation of 
the energy from the donor-acceptor system to the bath. 
The latter leads to the damping of coherent oscillations 
that would otherwise exist between the quantum states 
of the related molecule and, beyond a characteristic cou- 
pling strength, to incoherent dynamics of the electron 
transfer process. These considerations lead us to a dis- 
sipative two-site Hubbard model, a minimal model that 
captures the essential physics comprising correlations be- 
tween electrons and their coupling to the dissipative en- 
vironment. It is discussed in detail in Sec. [Til For a 
comparison to experimental results it has to be supple- 
mented by ab initio calculations of the parameters. 

The equilibrium properties of the model have been 
previously studied^ using the numerical renormaliza- 
tion group (NRG), and the real-time dynamics has been 
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using a Monte-Carlo technique at high 



temperatures where only incoherent transfer is present. 
In these Monte-Carlo calculations, the effective Coulomb 
interaction was chosen to be J7 e fF > and no electron- 
pair transfer has been reported. Two-electron transfer 
in a classical bath has been discussed in Ref. [2(| in the 
framework of three parabolic potential surfaces (for the 
four states D 2 ~A, D~A~ and DA 2 ~) as a function of 
a single reaction coordinate. A generalization to donor- 
bridge acceptor systems is given in Ref. (§l3l|. 

Although two-electron transfer was observed in some 
regimes of system parameters in the high-temperature 



limit, considering a classical bath, it seems reasonable to 
expect that, at least between identical centers, electron- 
pair tunneling processes are particularly important at 
temperatures corresponding to energies smaller than the 
effective energy difference between initial and interme- 
diate states U e g. At these temperatures single-electron 
transfer cannot be activated (see Section IVI[) . There- 
fore, we focus on the low-temperature regime where the 
transfer is dominated by nuclear tunneling and where 
the bosonic bath has to be treated quantum mechani- 
cally. Due to the nuclear tunneling the electron transfer 
rate is constant over a wide temperature range from zero 
temperature up to temperatures where thermal activa- 
tion becomes more important 32 . In this low-temperature 
regime, we employ the time-dependent NRG^ii 3 . (TD- 
NRG) which covers the whole parameter space from weak 
to strong dissipation. The NRG is an accurate approach 
to calculate thermodynamics and dynamical properties of 
quantum impurity modelsj 35 i 36 ' 37 i 38 For further details of 
the NRG we refer to the recent review 3 -^ on this method. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce the model. Its high temperature behavior ob- 
tained from the Marcus theory is described in Sec. IIIII 
Section HVl introduces the NRG method, its extension to 
non-equilibrium and its application to the present prob- 
lem. In order to gain a better understanding of the non- 
equilibrium dynamics presented in Sec. lVIi we summarize 
the equilibrium properties of the model in the Sec. El Wc 
present a detailed discussion of the real-time dynamics in 
Sec. IVII Therein, we focus on the time evolution of oc- 
cupation probabilities of the different electronic states as 
the key observables. In particular, when the dynamics 
can be described in terms of rate processes, the depen- 
dence of the single and electron-pair rate on the Coulomb 
repulsion parameters is analyzed. A summary of our re- 
sults is given in Sec. IVIII 



II. MODEL 

We consider a model of a two-electron/two-site system 
coupled to a bosonic bath. It is defined by the Hamilto- 
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where Ci a and c| CT denote annihilation and creation op- 
erators for fermions with spin er on the donor (i — D) and 
acceptor (i = A) sites. The Hamiltonian H e \ corresponds 
to an extended two-site Hubbard model, with on-site en- 
ergies Si, hopping matrix element A, on-site Coulomb 
repulsion U and an inter-site Coulomb repulsion V be- 
tween one electron on the donor and one electron on the 
acceptor. The difference U = U — V measures the ex- 
cess energy needed to move an electron between the two 
sites. Such a two-site Hubbard model without coupling 
to a bosonic bath has been investigated in the context of 
electron transfer in Ref. HjJ. 

The Hamiltonian Hb models the free bosonic bath, 
with boson creation and annihilation operators b n and b n , 
respectively. The electron-boson coupling term, H conp i, 
has the standard polaron form with the coupling con- 
stant for donor and acceptor given by gr>^n and gA^n, 
respectively. In what follows we set cd = — (a — § and 
gA = —gv = 1- The latter choice implies that the polar 
bath is coupled to the change in the electronic density 

Eo-( C L C Aa - C^CDa): 
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This two-site electron-boson Hamiltonian conserves the 
number of electrons J^ia c la c i>? an d tne sc i uare of the total 
spin S 2 as well as its z-component S z . The Hilbert space 
can therefore be divided into different subspaces. In the 
subspace with one electron and S z = 1/2, the model is 
equivalent to the spin-boson models. Here, we consider 
the subspace with two electrons and S z = which is 
spanned by the states |1) = | |j,0), |2) = | 1,T>, |3> = | T 
, J.}, and |4) = |0, Tl) with the notation \A, D) describing 
the occupation at the donor (D) and acceptor (A) sites. 
The four-dimensional basis in the two-electron subspace 
is displayed in Fig. [TJ We define the following observables 



do = |1)(1| 
d A = |4)(4| 
flAB = 1 2) (2 1 
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which measure the doubly occupancy d]j ((1a) on the 
donor (acceptor) site and tida the combined population 
of the states | f, [) and | |, t)- Note that in some works 
the states | tl>0) and 1 , J, T ) are referred to as localized 
states^. We call them doubly occupied states while the 
term localization is used below for the self-trapping mech- 
anism. 

Consider the 4x4 Hamiltonian matrix in the electronic 
subspace (M)ij = (i\H\j) — 1,...,4). Introducing 
the notation 
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FIG. 1: The four states of model (eq. (1)) for the symmet- 
ric case (e = 0). The energy difference between the doubly 
occupied donor (D 2 ~ A) or acceptor (DA 2 ~) and singly oc- 
cupied donor acceptor pair (D~ A~) depends on the effective 
renormalized interaction U c ff defined in eq. [5] 
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with U = U — V. Therefore, the dynamics of the sys- 
tem is governed by the energy difference U which re- 
places the on-site Coulomb repulsion U . If screening of 
the local Coulomb repulsion £ / 24 i 25 is sufficiently large, 
U changes its sign and become effectively attractive. A 
large inter-site Coulomb repulsion V favors an inhomo- 
geneous charge distribution. 

It is convenient to rewrite the diagonal matrix elements 
of the doubly occupied states in the form 

= c + U cK + Y,^n (bl + —) (bn + ^) , (6) 



and 

(4|ff|4> =-e + U cS |V,„ (bl ^) (bn - ^) .(7) 

Compared with the matrix elements of states correspond- 
ing to D~A~ 



{2\H\2} = {3\H\3} =J2"nb n b n , 



(8) 



we can easily see that the electron-boson coupling gener- 
ates an effective renormalized interaction 



(9) 



The renormalized interaction C/ C ff determines the energy 
difference between D 2 ~ A (DA 2 ~) and D~A~, and con- 
stitutes the only Coulomb interaction parameter in the 
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present model. The renormalization stems from a boson- 
induced effective electron-electron interaction, already fa- 
miliar from the Holstein model 2 ^. Note that an artifi- 
cial energy shift is present in the single-electron subspace 
(spin-boson model)£, however the two states | |) and | j) 
arc shifted in the same direction which can be handled 
by resetting the zero of energy. 

In analogy to the spin-boson mode&£, the coupling of 
the electrons to the bath degrees of freedom is completely 
specified by the bath spectral function 
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The spectral function characterizes the bath and the 
system-bath coupling, and can be related to the classical 
reorganization energy^ (classical in terms of boson de- 
grees of freedom) which measures the energy relaxation 
that follows a sudden electronic transition. The one- 
electron transfer and the correlated two-electron trans- 
fer are associated with reorganization energies E a \ and 
E a2 , respectively. For a single-electron transfer, e.g., 
D 2 ~ A —¥ D~A~ the reorganization energy E a \ is given 
by£ 



A 2 

E a l = > — = 

^ U)„. 



duj J(lu) 
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and the corresponding energy for a correlated two- 
electron transfer (D 2 ~A — » DA 2 ~) is E a2 = 4:E al . 

The model we are considering here is completely spec- 
ified by the parameters A, a, U, e and the bosonic spec- 
tral function. In the molecular electron transfer problem 
the latter function reflects intramolecular vibrations and 
the solvent (e.g., water or protein) or environment. Its 
solvent component can be estimated from the solvent di- 
electric properties or a classical molecular dynamics sim- 
ulation. In the present paper we assume an Ohmic bath 
model: 



FIG. 2: Potential surfaces for the different states of the model 
in the Marcus theory for U e g > 0, e = 0. The mimima of the 
states | |,|) and | J.,|) (D~A~) are set to the origin while 
those parabolas that correspond to the doubly occupied states 
| T|,0) (D 2 ~A) and |0, Ti) (DA 2 ~) are shifted. Note that in 
the case displayed here the transfer D 2 ~ A — > D~ A" is in the 
"inverted regime" . 



non-adiabatic limit) are also available^, however for sim- 
plicity we limit ourselves in what follows to the classical 
Marcus description. The Marcus electron transfer rate 
can be evaluated for any amount of transferred charge: 
the latter just determines the renormalized potential sur- 
face parameters that enter the rate expression. Single- 
electron transition rates are given by 
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with a cut-off at energy uj c . This choice yields the re- 
organization energy E a \ = 2auo c and the energy shift 
U c ff = U — 2au> c . All parameters and physical quanti- 
ties are defined in units of lo c . Its order of magnitude for 
the intermolecular mode spectrum of a polar solvent is 
O.leV. 



\D~A-^DA 2 
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In the case A <C \ U c s | second-order processes are possible 
that involve only virtual occupations of the states D~A~ 
leading to rates for an electron pair. 



i pair 
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III. THE HIGH-TEMPERATURE LIMIT: 
MARCUS THEORY 

In the high-temperature, limit electron transfer is usu- 
ally described using Marcus theory^ as a rate process 
within classical transition state theory. Extensions that 
take into account the quantum nature of the nuclear mo- 
tion in the weak electronic coupling limit (the so called 



A 4 (2e + E a2 ) 

K [DA 2 ~^D 2 -A] f T 2 6 ■ 

The interplay between sequential and concerted two- 
electron transfer (in the limit of a classical bath with 
a single mode or a single reaction coordinate) can be 
seen from these expressions. In the following we re- 
strict ourselves to the symmetric case (e — 0). Starting 
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with the initial state D 2 A, we expect concerted two- 
electron transfer in the Marcus regime when the rate 
k [D*-A->DA 2 -] is l ar g er tl:l an the rate k s ^2- A ^ D - A -] of 
the first step of the sequential process which is the case 
when \U cS \ > T and \U cS \ > E al as well as E aX < T. 

In a parameter region where sequential transfer dom- 
inates the rates k \"jj2— a >d~ a— ] and fcr^— ^4— ^z^a 2- ] a ^ 

well as the corresponding backward rates show a non- 
monotonic behavior and an inverted regime dependent 
on the effective Coulomb interaction U eB (see Fig. [2]). 

For incoherent transfer processes (which may happen 
at large temperatures and for a strong coupling to the 
bosonic bath), a description of the population dynamics 
by kinetic equations determined by the rates is given by 

d D (t) = - (^ n £° A ^ D - A -] + h^- A _^ DAa -yj d D {t) 

+ ^ S [D-A-^D 2 -A] TlDA ( t ) + k [DA 2 -^D 2 -A] dA ( t ) ' 
h DA (t) = - (fc[^ S A-^,D^-] + k [D-A-^D*-A]) n DA(t) 
+ 2fc[ D 2°- A-^D- + ^[DA 2 - ^D- A-]^ A ^) ' 
d A(t) = - {k S [DA^-^D-A-] + k "[DA 2 -~>D 2 -A]) ^ W 

+k [D-A~ -,DA 2 -] nDA (~ t ) + ^"-A^M^l^W' 

(19) 

where do and g?a are the probabilities to have two elec- 
trons on the donor and acceptor, respectively, uda is 
the combined population of the states | T, j) and | j, |). 
For the initial condition d[)(t = 0) = 1 we obtain 
Tl |T-l) = "li-T)- F° r the unbiased Hamiltonian (e = 0), 

i- single — jingle m J ^single 

^[DA 2 -^D-A-] — ^[D 2 - A^D- A~] dI1U K [D~ A- -^DA 2 ~] ~ 

k [D S A-^D 2 -A] nlUSt hold - 

These kinetic equations can be solved in the high- 
temperature regime using the Marcus rates from above. 
In Sec lVIl we have used these equations to extract 
the low-temperature transition rates by fitting the non- 
equilibrium dynamics of du(t),dA(t) and nDA(t) calcu- 
lated in the incoherent regime with the time-dependent 
NRG. 

For t — > 00 the equilibrium states (dA)eq, (rfo)eq and 
(^DA)eq are reached, where (djj) eq — (e^cq- It follows 

that , {dD) ^ = k[ °~it~^ DA2 ~ ] which is according to the 

[DA 2 - —*D~A~] 

Marcus rates .^-"^ — e^ cff / T . Therefore, in the classi- 
cal limit we arrive at 

(d D )f a = -° ,S ■ (20) 

N /C< 1 e U/T _|_ j y ' 

With the help of the kinetic equations we can de- 
scribe concerted two-electron transfer, a purely sequen- 
tial single-electron transfer as well as a combined pro- 
cess which shows first a pair transfer which is fol- 
lowed by a single-electron transfer. As long as the 
single-electron transfer rates are small (fcfj^fj a—>d- a-] < ~ 



k [D 2 -A^DA 2 -]) and ( d D) cq = (eU)cq ~ 0.5 the state 
D~ A~ is only weakly populated and noA(t) is con- 
stant and close to zero. The dynamics is dominated by 
an electron pair transfer. The combined process is ex- 
pected if (dr))eq = (cU)cq < 0.5. First the population 
dA rises quickly while n^A stays close to zero. Later a 
slow increase of uda to its equilibrium is observed. For 

k fp*-A->D-A-] > k [D 2 -A^DA 2 -] tne transfer is purely 
sequential. 



IV. THE LOW-TEMPERATURE LIMIT: THE 
NUMERICAL RENORMALIZATION GROUP 

At low temperature, the quantum generalization of the 
Marcus theory replaces the classical environment by a 
bath of non-interacting bosonic degrees of freedom. Very 
early on, the " non-adiabatic" weak coupling limit was 
investigated^ The strong coupling limit of such a model 
has been adressed using the non-interacting blib approx- 
imation (NIBA)£, path integral methods^ and recently 
also by the numerical renormalization group which we 
employ in this paper. 

Originally the NRG was invented by Wilson for a 
fermionic bath to solve the Kondo problemj 35 ' 36 The 
fcrmionic NRG is a standard and very powerful tool to 
investigate complex quantum impurity problems^ The 
method was recently extended to treat quantum im- 
purities coupled to a bosonic bat bj L S S., to a combina- 
tion of fermionic and bosonic baths 4 ^, and to the calcu- 
lation of real-time dynamics out of equilibrium i 33 i 34 i 43 
The non-perturbative NRG approach has been suc- 
cessfully applied to arbitrary electron-bath coupling 
strengthsffl^^i 



A. Equilibrium NRG 

The numerical renormalization group achieves the sep- 
aration of energy scales by logarithmic discretization of 
the energy continuum into intervals [A - ( m+1 'ui c , A~ m ui c ], 
m G Ao, defining the discretization parameter A > 1. 
Only one single mode of each interval couples directly 
to the quantum impurity, indicated by the circles in 
Fig. GH a ) . This discrete representation of the continuum 
is mapped onto a semi-infinite tight-binding chain us- 
ing an exact unitary transformation. Hereby, the quan- 
tum impurity couples only to the very first chain site as 
depicted in Fig. (3Jb). The tight-binding parameters t n 
linking consecutive sites of the chain m and m + 1 fall 
off exponentially as t m ~ A~ m . Each bosonic chain site 
is viewed as representative of an energy shell since its 
energy w m also decreases as w m ~ A~ m establishing an 
energy hierarchy. Both ensures that mode coupling can 
only occur between neighboring energy shells which is es- 
sential for the application of the renormalization group 
procedure. To this end, the renormalization group trans- 
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FIG. 3: Scheme of the bosonic NRG. (a) The bosonic energy 
continuum is discretized on a logarithmic mesh using a pa- 
rameter A > 1. Only a single bosonic mode in each interval 
[A~ < - n+1 ' > uj c , A _n tj c ] - visualized by the circles - couples di- 
rectly to the electronic subsystem, (b) This discretized model 
is mapped exactly onto a tight-binding chain via a unitary 
transformation^^ only the first chain site couples directly 
to the donor-acceptor system. The hopping t„ between neigh- 
boring bosonic sites decreases exponentially with the distance 
from the donor- accept or system, i.e. t n oc A~ n . The energy 
spectrum of the Hamiltonian is calculated by successively ap- 
plying the renormalization group (RG) transformation (|21|l . 
diagonalizing the new Hamiltonian and rescaling the spec- 
trum as depicted schematically in panel (c) for the sequence 
of Hamiltonians H m to H m +2,- After each iteration only the 
N s eigenstates of site m + 1 with the lowest energies are kept. 
This truncation is depicted by a horizontal dashed line. 



Hamiltonians under the RG transformation R[H m ]. 

The RG transformation (f2"Tj) is used to set up and iter- 
atively diagonalize the sequence of Hamiltonians H n . In 
the first step, only the electronic donor-acceptor system 
coupling to the single bosonic site to = is considered. 
It turns out to be sufficien t 37 ' 38 ! 39 to include only the iVf, 
lowest lying bosonic states, where Nb takes typical values 
of 8 — 12. The reason for that is quite subtle: the cou- 
pling between different sites decays exponentially and is 
restricted to nearest-neighbor coupling by construction, 
both essential for the RG procedure. In each successive 
step (i) a finite number of Nb bosonic states of the next 
site to+1 are added, (ii) the Hamiltonian matrices are di- 
agonalized and (iii) only the lowest N s states are retained 
in each iteration. The discarding of high-energy states 
is justified by the Boltzmannian form of the equilibrium 
density operator when simultaneously the temperature is 
lowered in each iteration step to the order T m tx A~ m w c . 

To illustrate the procedure, the lowest-lying energies of 
the Hamiltonian H m to H m+ ^ are schematically depicted 
in panel (c) of Fig. [3] We typically use Nb > 8 and 
keep about N s = 100 states after each iteration using a 
discretization parameter A = 2. 

Denoting the set of low- lying eigenstates by |r)jv and 
the corresponding eigenvalues E r (N) tx 0(1) at iteration 
N, the equilibrium density matrix po is given^ 9 - by 



Po 



1 

Zn 



e-~ 0E ?\r) NN {r\ , 



(22) 



where Zn = J2 r e ,JKr an< ^ P are °f the order O(l), 
such that Tjy = w c A~ N / (3. The thermodynamic expec- 
tation value of each local observable O is accessible at 
each temperature T/v by the trace 



(6>oq = Tr 



PoO 



1 

Zn 



N (r\6\r) N .(23) 



The procedure described above turns out to be very ac- 
curate because the couplings t m between the bosonic 
sites along the chain are falling off exponentially, so 
that the rest of the semi-infinite chain contributes only 
perturbativel y 35 ' 39 at each iteration to while contribu- 
tions from the discarded high-energy states are exponen- 
tially suppressed by the Boltzmann factor. 



formation R[H] reads 

H m +i — R[H m ] = Ai? TO + A' Tl+1 (i m flj„ffl m +l 

+ t m a\n +1 a m + Wm,ain+l a m+l) ' ( 21 ) 

where H m is the Hamiltonian of a finite chain up to 
the site to - as depicted in Fig. E^b). The annihilation 
(creation) operators of site m are denoted by a m (a^) 
and w m labels the energy of the bosonic mode of site to. 
Note that the rescaling of the Hamiltonian H m by A en- 
sures the invariance of the energy spectrum of fixed point 



B. Time dependent NRG 

While the equilibrium properties are fully determined 
by the energy spectrum of the Hamiltonian, the non- 
equilibrium dynamics requires two conditions: the ini- 
tial condition encoded in the many-body density op- 
erator po and the Hamiltonian H f which governs its 
time-evolution. For a time-independent Hamiltonian, 
the density operator evolves according to p(t > 0) = 
e -iH tp Qe tH *. All time-dependent expectation values 
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(0)(t) are given by 

(6}(t) = Tr\p(t)6 



Tr 



-iHH iHH 



p e l " £ 



(24) 



We obtain the density operator po from an independent 
NRG run using a suitable initial Hamiltonian H 1 . By 
choosing appropriate parameters in H 1 , we prepare the 
system such that (for the calculations presented in this 
paper) the two electrons are located on the donor site 
and the acceptor site is empty. 

In general, the initial density operator po contains 
states which are most likely superpositions of excited 
states of H l . For the calculation of the real-time dynam- 
ics of electron-transfer reactions it is therefore not suffi- 
cient to take into account only the retained states of the 
Hamiltonian H f obtained from an NRG procedure. The 
recently developed time-dependent NRG (TD-NRG) 3 ^ 
circumvents this problem by including contributions from 
all states. It turns out that the set of all discarded states 
eliminated during the NRG procedure form a complete 
basis se t 33 ' 34 of the Wilson chain which is also an approx- 
imate eigenbasis of the Hamiltonian. Using this complete 
basis, it was show n 33 ' 34 that eq. transforms into the 
central equation of the TD-NRG for the temperature Tjv 

N trim 

< 6 >W = E E e^-^O^pl^m) , (25) 



where 0" l s — (r;m\0\s;m) are the matrix elements of 

any operator O of the electronic subsystem at iteration 
m, and E™,E™ are the eigenenergies of the eigenstates 
\r;m) and \s;m) of H l m . At each iteration m, the chain 
is formally partitioned into a "system" part on which 
the Hamiltonian H m acts exclusively and an environment 
part formed by the bosonic sites m + 1 to N. Tracing 
out these environmental degrees of freedom e yields the 
reduced density matri x 33 ' 34 



pl c ,r ( m ) = E( s ' e ; m \Po\ r i e ; m ) 



(26) 



at iteration m, where po is given by (f2"2"| using H 1 . The 
restricted sum X^'s" m ec f P^)) implies that at least one 
of the states r and s is discarded at iteration m. Exci- 
tations involving only kept states contribute at a later 
iteration and must be excluded from the sum. 

As a consequence, all energy shells m contribute to the 
time evolution: the short time dynamics is governed by 
the high energy states while the long time behavior is 
determined by the low lying excitations. Dephasing and 
dissipation is encoded in the phase factors ~ E » as 
well as the reduced density matrix p T f~(m). 

Discrctizing the bath continuum will lead to finite-size 
oscillations of the real-time dynamics around the contin- 
uum solution and deviations of expectation values from 
the true equilibrium at long time scales. In order to sep- 
arate the unphysical finite-size oscillations from the true 



continuum behavior, we average over different bath dis- 
cretization schemes using Oliveira's z-averaging (for de- 
tails see Refs. MM)- 

We average over 8 different bath 
discretizations in our calculation. 



V. EQUILIBRIUM PROPERTIES 

In order to gain a better understanding of the non- 
equilibrium dynamics presented in Sec. IVI1 we briefly 
summarize the equilibrium properties of the model given 
by Eq. (JT|). It has been analyzed in Ref. (28|, where 
self-trapping (localization) in the single and two-electron 
subspace was found. 

We start with the phase diagram of the two-site model, 
as shown in Fig. 0J Only for e = a quantum phase tran- 
sition of Kosterlitz-Thouless type separates a localized 
phase for a > a c from a delocalized phase for a < a c . 
We plot the phase boundaries between localized and de- 
localized phases in the a-[/-plane, both for single- and 
two-electron subspaces (grey and black line in Fig. [U re- 
spectively). 

For the single-electron subspace, the Coulomb repul- 
sion is irrelevant, and the phase boundary does not de- 
pend on U. The value of the critical coupling strength, 
a c , is identical to those of the corresponding spin-boson 
model. The critical value^^ of a c depends on the tun- 
neling rate A and reaches a c — 1 for A — » 0. 

The phase boundary for the two-electron subspace does 
depend on U, which has drastic consequences for the elec- 
tron transfer process. Imagine that, by a suitable choice 
of parameters, the system is placed between the two 
phase boundaries above the single-electron (grey line) 
and below the two-electron phase boundary (black line) 
in the area indicated by I in Fig. 3J Then the system 
would be in the localized phase in the single-electron sub- 
space. However, one additional second electron immedi- 
ately places the system in the delocalized phase, and one 
or even both electrons can be transfered. Similarly, a 
second electron added to the system in the parameter 
regime of area II shows the opposite behavior: both elec- 
trons get localized although a single electron could be 
transfered. 

Note the different values of the a c 's even for U — 
in the single and the two-electron subspace: the cou- 
pling of the donor/acceptor system to the bath induces 
an effective attractive Coulomb interaction U c s = —2auj c 
between the electrons. On the localized side of the 
transition, the electron tunneling A is renormalized to 
zero, so that an electron transfer is clearly absent in this 
regime. This statement holds only for Ohmic dissipation, 
on which we focus here; deep in the sub-ohmic regime, 
coherent oscillations have been recently observed even in 
the localized phase, see Ref. 0- 

Figure [5] shows results for the double occupation prob- 
ability as a function of the electron-bath coupling a for 
different U calculated with the equilibrium NRG. For 
the symmetric model considered here, the equilibrium 
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FIG. 4: Zero-temperature phase diagram of the model eq. ((TJ 
for e = and A = 0.1oj c . The critical dissipation strength 
a c is plotted as a function of (7 in the two-particle subspace 
(black line) and in the single-particle subspace (grey line), 
respectively. 



0.5 



0.4 



0.3 



A 
T3 

V 0.2 



0.1 



0.0 



h-ri — ' — ' u — *- 
0=-a> / 

■ 


•— ^x__jy — I 1 1 — 


. " " 

/ 


: / 

* u=o 




■ 

/ : 




0=co c / 

■ 

/ 




a ■ •■ ■ ' 

i . i 








0.0 0.1 0.2 0.3 0.4 0.5 0.6 

a 



FIG. 5: Low-temperature equilibrium probability for double 
occupancy of donor and acceptor (d) eq for A = 0.1u c , e = 0, 
as a function of a for U = —ui c , and lo c . In the limit of a = 
the dependence of (d) oq on A and U is given analytically in 
eq. (27J). 



probabilities for the double occupation on donor and 



acceptor sites are equal: (d\} cq = (du}e 



ing the observables defined in Eq. ([3]). The probabil- 
ity of having two electrons at different sites is given by 

("DA)cq = 1 - 2(d) C q. 

The average double occupancy (d) eq decreases with 
increasing effective Coulomb repulsion U and increases 



FIG. 6: Equilibrium probability for double occupancy of 
donor and acceptor (d) eq as a function of temperature T for 
A = O.lWc, Ccff = O.luc, a = 0.04 (circles, dashed line) and 
or A = O.Iluc, U c ff = U) c , a — 0.04 (squares, solid line). For 
comparison the "high temperature" result eq. (|20[) is shown 
for (7 e fr = 0.1aj c (dashed line) and U c s = lo c (solid line). 



with increasing a. This can understood in terms of the 
effective Coulomb interaction U e g — U — 2aw c , renormal- 
ized due to the coupling to the bosonic bath. 

The delocalization/localization phase transition occurs 
when (d)cq — > 0.5, as can be seen by comparing Fig.[4]and 
Fig. [5] For U e s < and U e s >• A, we are able to project 
out the D~A~ excited states. Then our model maps 
on a spin-boson model with an effective hopping A/U^ s 
between the states D 2 ~A and DA 2 ~ . The dynamics will 
be governed by electron pairs if D 2 ~ A or DA 2 ~ are the 
initial states. 

The double occupancy {d) cq is calculated analytically 
for a — and arbitrary A and U. For T — > 0, the (d) eq 
approaches 



(d) 



4A 2 



eq 



VU 2 + 16A 2 (U + y/U 2 + 16A 2 ) 



(27) 



while in the opposite limit, T — * oo, we obtain (d) cq — ► 
0.25. The low-temperature limit (|27|) is included as end- 
points of the curves in Fig. [5] 

Let us now turn to the temperature dependence of 
(d) q. Figure [5] shows results for temperatures between 
T = 0.004w c and T = 0.2u> c for several choices of 
model parameters. Our calculations imply an indepen- 
dent check of the correct t — > oo behavior in the next 
section. Additionally, we can make connection to the 
high-temperature results of section IIIII For tempera- 
tures T <ti U c g the double occupancy (d) Qq is constant 
as expected from quantum statistics but deviates dras- 
tically from the predictions of the Marcus theory given 
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by Eq. ([20]) . The double occupancy (djj) cq calculated 
with the NRG approaches the value 0.5/(1 + e u " !! ^ T ) for 
U c fi ~ T. This result indicates that for U c g > T Marcus 
theory is not applicable while low temperature methods 
like the NRG are valid. 



VI. NON-EQUILIBRIUM DYNAMICS 

We employ the time-dependent NRG to evaluate the 
low-temperature time evolution of the local occupancies 
using Eq. (f25|) and investigate the influence of differ- 
ent Coulomb interactions U, single-electron hopping ma- 
trix elements A, couplings between the electronic sys- 
tem to the bosonic bath a and temperatures T between 
T = 3 • 10~ 8 w c and T = 0.125w c . The donor/acceptor 
sub-system is initially prepared in a state with the two 
electrons placed on the donor site and evolves according 
to Hamiltonian |T|). We calculate the time-dependent 
expectation values dr>(t) = (do)(t), dx(t) = ((1a) (t) and 
nuA(t) — (fiDA}(t) using Eq. These expectation 

values are related at any time by the completeness rela- 
tion d,D{t) + d,A(t) + nDA(t) = 1. The time evolution of 
riDA(t) serves as criterion to distinguish between direct 
two-electron transfer and two consecutive one-electron 
steps. If n D A{t) remains close to zero or stays constant 
throughout the electron transfer process, the two states 
D~ A~ are only virtually occupied, and concerted two- 
electron transfer is observed. A significant increase of 
nDA(t) as function of time is taken as an indication of 
step-by-step single-electron transfer. 

In the absence of the electron-boson coupling (a = 0), 
the dynamics is fully determined by the dynamics of the 
four eigenstates of H e i. In the limit \U\ 3> A we obtain 



dD-.At) 
n DA (t) ■■ 



1 2A 2 2A 2 

2 ~ ~U 2 + Jp 
4A 2 4A 2 



cos 







)±icos 


;f<) 







u 2 u 2 



■ cos 



(28) 



while in the limit of U = 
3 1 
8 + 8 



ID-A 



(t) 



cos(4At) ± - cos(2At) , 



2 2 

n D A{t) = -- -cos(4At). 

o o 



(29) 



A finite value of the coupling, a ^ 0, gives rise to 
damping of those coherent oscillations. Furthermore, the 
Coulomb interaction is renormalized to U c s = U — 2auj c . 
For U c s < 0, the states D 2 ~ A and DA 2 ~ are energet- 
ically favored. The two intermediate states D~A~ are 
only virtually occupied for \U e g\ 3> A, T, similar to the 
superexchange process. 1 This regime can be described by 
a spin-boson model with an effective interstate coupling 
A e ff « 4A 2 /f/ e fj. The spin-boson model has three dy- 
namical regimes.—. 

For a smaller than some characteristic value it exhibits 
damped coherent oscillations between the two states. If 
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FIG. 7: Low-temperature population probabilities P(t) = 
do(t) (thick black line), d,A(t) (thin black line) and riDA(t) 
(gray line) as functions of time. The parameters are U = — lj c , 
A = O.lWc, e = and T = 3 ■ 10~ 8 lu c . The coupling a in- 
creases from the upper panel a = (f7 e ff = — w c ), the middle 
panel a — 0.04 (f7 e ff = — 1.08o; c ) to the lower panel a = 0.16 
(U eH = -1.32u> c ). 



a is larger than this value the oscillations disappear and 
the kinetics is dominated by a relaxation process. Here, 
rates can be defined and the population probabilities can 
be fitted with the kinetic equations (TT5)) . For a further 
increase of a, the electronic system shows the onset of 
localization (for T — > 0) and does not evolve towards the 



10 




50 100 150 200 ' 50 100 150 200 



FIG. 8: Low-temperature population probabilities P(t) — 
do(t) (thick black line), d,A(t) (thin black line) and noA(t) 
(gray line) as functions of time. The parameters are U = lo c , 
A = 0.1w c , e = and T = 3 • 10" 8 lj c . The coupling to 
the bosonic bath increases from panel (a) a — 0(t/ c ff = We)j 
panel (b) a = 0.02 (U eS = 0.96w c ), panel (c) a = 0.52 
(U c ff = -0.04w c ) to panel (d) a = 0.55 (U e g = -0.kj c ). 

other (acceptor) site. 

In Fig. [7|we plot the low-temperature population prob- 
abilities do(t), d.A(i) and noAit) for U = —uj c , A = 
Q.lw c , T = 3 ■ 10~ 8 w c and different couplings a. For 
a = (upper panel) the oscillations have two frequencies 
(see eq. (|28p). The electron pair oscillates from donor 
to acceptor with the small frequency 4A 2 /{7 whereas the 
fast oscillations with frequency U characterize the virtual 
occupation of the high lying states (D~A~). An increase 
of a leads to damping of the oscillations (middle panel) 
and relaxation (lower panel). At about a = 0.3 the elec- 
tron pair gets self-trapped and the system shows a phase 
transition to the localized phase at T=0 (see Fig.[5J). The 
configuration D~ A~ is seen not to be involved in the dy- 
namics as no a is very small and without ascending slope. 
Since A <C U e s the state D~A~ cannot be populated as 
long U cS > T. 

A more complicated behavior is expected within the 
four accessible electronic states when U e g 3> A > 0. In 
this case the delocalized states D~A~ have the lowest 
energy, and sequential transfer is required to reach the 
equilibrium state. Pair transfer occurs on a smaller time 
scale. Thus, a combined pair and sequential transfer on 
two different time scales governs the dynamics for these 
parameters. 

The four panels in Fig. [5] depict the time evolution 
of the occupation probabilities do{t), dA(t) and noA{i) 
for U — u) c , A = 0.1o; c and four different values of a: 




FIG. 9: Low temperature population probabilities P(t) = 
doit) (black line) and riDA(t) (gray line) as functions of time 
for a = 0.04 (full line), a = 0.36 (dashed line). The effective 
energy difference between the states D 2 ~ A and D~ A~ is kept 
constant U c a = ui c . The other parameters are A = 0.1ui c and 

£ = 0. 
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FIG. 10: Low temperature population probabilities P(t) = 
do(t) (thick black line), dA(t) (thin black line) and noA(t) 
(gray line) as functions of time. The parameters are U = 
0, A = 0.1o; c , £ = and T = 3 • 10" 8 lj c . (a) a = 0((/ cff = 0), 
(b) a = 0.01((7 cff = -0.02o; c ), (c) a = 0.04((7 cff = -0.08w c ) 
and (d) a = 0.1((7 cff = -0.2w e ). 



a = 0, 0.02, 0.52, 0.54. The undamped coherent oscilla- 
tions of panel (a) decay exponentially for small damping 
depicted in panel (b). Increasing a further yields a fi- 
nite population of the states D~ A~ : sequential transfer 
becomes the main process, as shown in panels (c) and 
(d). The crossover from a combined pair-transfer and 
slow single-electron transfer (panel(b)) to purely sequen- 
tial transfer (panel (c) and panel (d)) with a complex 
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dynamics is due to a combined effect of dissipation and 
decrease of the effective energy difference between the 
donor/acceptor states. An even larger a leads to a neg- 
ative U e s and a very slow transfer until the onset of lo- 
calization at a c which is not shown here. 

To separate the influence of dissipation from the renor- 
malization of U due to the coupling to the bosonic 
bath we plot doit) and Uda{P) for a constant effective 
U e a = oj c and different coupling a in Fig.0 The dynam- 
ics changes from pair transfer with a slow increase of the 
single occupancy at a = 0.04 (due to the low-lying states 
D~A~) to incoherent relaxation and sequential transfer 
for a — 0.36. As long as E a i = 2au c <C U c s, pair trans- 
fer is observed on a short time-scale. For E a \ > U c g only 
one electron is transferred and the system relaxes rapidly 
into its equilibrium state D~ A~~ without any short time 
concerted pair transfer. 

In Fig. [10] the evolution of the dynamics is shown for 
U = and increasing a. The doubly occupied states 
are the ground states of the donor/acceptor system for 
finite a since U c ff = —2auj c < 0. With increasing a, the 
amplitude of coherent oscillations acquire a small damp- 
ing. In addition, pair transfer is favored and nj)A(t) de- 
creases. The simple damped oscillations are replaced by 
a much more complex dynamics comprising of strongly 
renormalized oscillation frequency and a strong damping 
for a = 0.04. At about a — 0.36 - not shown here - the 
critical coupling a c is reached and the system is localized. 

Next we study the effect of changing U at constant sys- 
tem bath coupling a = 0.04 (Fig. [TTJ), a = 0.16 (Fig.[T2|) 
and a = 0.36 (Fig.rjJ) and A = 0.1w c . 

In the lower damping case (Fig. [TTj) , the transfer is re- 
flected by damped electron pair oscillations for U = —ui c 
in Fig. [TTJ (a). Increasing U — — 0.5w c in Fig. [TTJ (b) leads 
to an increase of the population probability of D~ A" and 
to a change of the fast oscillations whit an approximate 
frequency of U c s. When U becomes positive U = 0.5w c 
(Fig. [TTJ (c)) the single-electron transfer becomes fast and 
the main process unless U D s is not too large. In fact at 
U = lu c the rate from D 2 ~ A to D~A~ becomes smaller 
(Fig. [TTJ (d)) and additional electron pair transfer is ob- 
served. The graphs Figs. [TTJ a) and (d) can be under- 
stood in terms of Eq. {28j since \ A/U\ <C 1 and a = 0.04 
is small. By the weak coupling to the environment, A 
is slightly reduced, and the oscillation amplitude decays 
exponentially. The difference between the two panels (a) 
and (d) arises from (i) U c ff = U — 2auj c instead of the \U\ 
entering Eq. (|28|) and (ii) from the dissipation which fa- 
vors the relaxation into the new thermodynamic ground 
state: while the oscillation frequencies are roughly the 
same for \U\ = u> c , the delocalized states have a lower 
energy in Fig. ITTT d) so that riDA(t) has to increase to 
its new equilibrium value. The approximations made in 
Eq. (|28|) do not hold any longer for the parameters in 
Figs. [JjTb) and [TlTc). The electronic dynamics is gov- 
erned by additional frequencies and becomes more com- 
plex. However, the results can still be analyzed and un- 
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FIG. 11: Low temperature population probabilities P(t) = 
do(t) (thick black line), dA(t) (thin black line) and noA(t) 
(gray line) as functions of time. The parameters are a — 0.04 
and U = —u>c (panel (a)), U = — 0.5w c (panel (b)), U = 0.5w c 
(panel (c)) and U — u> c (panel (d)) 
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FIG. 12: Upper panel: Low-temperature population proba- 
bilities Pit) = doit) (thick black line), d,A(t) (thin black line) 
and n,DA(t) (gray line) as functions of time. U = — 0.9u> c 
and — 1.5ui c from bottom to top for do as well as from top 
to bottom for (La- Lower panel: Electron pair rate fc pa i r (for 
the transfer from D 2 ~ A — * DA 2 ~) as a function of U e s- The 
parameters for both panels are a — 0.16, T ~ 3 ■ 10 _8 oj c , 
A = O.Iluc and e = 0. Inset: Energy levels of states D 2 ~ A, 
D~A~ and DA 2 ~ . 



derstood within the analytical results of drj(t), d,A(t) and 
nDA(t) for a — ► 0. 

When the coupling a is increased to a — 0.16, a dif- 
ferent picture emerges. Very high frequency oscillations 
with a small amplitude are superimposed on a slowly de- 
caying do{t) depicted in in the upper panel Fig. 1121 Av- 
eraging over those oscillations, we can fit the population 
probabilities to the kinetic equations (jTUJ) . By this proce- 
dure, we extract the phenomenological rates as function 
of t/ e ff for fixed a = 0.16. As shown in the lower panel 
of Fig. [T2]the concerted transfer rate kTjj2- a-^da 2 ~] in- 
creases with increasing t/ e ff (U e g < 0). This was expected 
from the rate equation (JT7J) in the classical limit. 

The transfer is found to be incoherent and sequential in 
the higher damped case a = 0.36 for not too large U. The 
population probabilities are shown for U = 0.6w c , 1.7u> c 
and 2.5w c in the upper panel of Fig. [T3J By fitting the 
curves with the help of the kinetic equations, eq. (fT!))) . 
we obtain the rate of the single-electron transfer D 2 ~ A 
to D~ A~ which is a non-monotonic function of U e g with 
a maximum at U e g = E a \ « 0.72u; c (see lower panel). It 
is plotted together with the Marcus rate at T = 0.008w c 
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FIG. 13: Upper panel: Low-temperature population proba- 
bilities P(t) — do{t) (thick black line), dA(t) (thin black line) 
and riDAit) (gray line) as functions of time. For dA from 
top to bottom U — 0.6oj c , 1.7w c and 2.5lj c . The other pa- 
rameters are a = 0.36, T ~ 3 • 10~ 8 lo c , A = 0.1u; c and e = 0. 
Middle panel: Low-temperature population probabilities with 
U — 5uj c - Lower panel: Single-electron rate fc S m g ie for the 
transfer from D 2 ~ A — > D~ A~ deduced by fitting the popula- 
tion probabilities with the kinetic equations eq. (|19p (squares) 
and Marcus rate eq. (|13[) (full line) with T = 0.008o; c as a 
function of the on-site Coulomb repulsion. The Marcus rate is 
normalized so that both curves have the same maximal rate. 
Inset: Energy levels of states D 2 ~ A, D~ A~ and DA 2 ' . 



13 



T<0.02 co 




FIG. 14: Population probability uda of the state D~A~ as 
a function of time t for temperatures between T < 0.02oj c 
and T = 0.125oj c . The parameters are U = — 0.01cj c , A = 
O.OOlwc, a = 0.03. 

(For varying temperatures we found that the fitted rate 
is approximately constant for temperatures T < 0.008a-> c 
in the considered parameter space.). Although the quali- 
tative behavior is captured by the Marcus rate the asym- 
metric shape of the NRG result is more realistic in the 
nuclear tunneling regime. As U increases further the se- 
quential transfer becomes negligible in the inverted re- 
gion. As a matter of fact, an increasing value of U shifts 
the system away from the phase transition line deeper 
into the delocalized phase as can be seen in the equi- 
librium phase diagram Fig. 2] Here, the dynamics is 
dominated by coherent pair oscillations with a very small 
frequency displayed for U — 5u c in the middle panel of 
FigrU 

Finally, the effect of temperature is studied in Fig. [Til 
where U = — 0.01o> c , A = 0.001w c . The temperature is 
varied from 3 • 1(TV to 0.125w c . For T = 3 ■ lCTV 
to T < 0.02w c the population probability is tempera- 
ture independent. As long as U c g > T, pair transfer is 
observed (the probability of D~A~ stays constant). As 
T > U c ff the states D~A~ are seen to contribute and are 
thermally populated. 



VII. SUMMARY AND CONCLUSION 

In this paper, we have studied the electron transfer 
properties of two excess electrons in a redox system mod- 
eled as a dissipative two-site Hubbard model - a model 
which can be viewed as the simplest generalization of 
the spin-boson model to include many-particle effects. 
These many-particle effects are due to on-site and inter- 
site Coulomb interactions, U and V respectively, as well 
as the effective interactions induced by the coupling to a 



\U cft \ < A, 


single- 


electron transfer 


T> U cff 


single- 


electron transfer 


UeB > 0: 






U c g < E a l 


single- 


electron transfer 



k[D 2 - a~ >d-a~] fester and 
k [D - A ~^ D - A 2- ] slower 
with increasing [/ e ff 
Ucft > E a \ single-electron transfer 

k[D 2 -A^D-A~] 

and k [D - A -^ D - A 2-] 
slower with increasing U e s 

U cB > A, \Ueft\ > E al , 

T < | [/eff | electron-pair transfer 

(in addition slow 
single-electron transfer) 

UeB < 0: 

\U e s\ < E a i single-electron transfer 

fc[ D - A -_ D - A 2-j faster and 
k [D 2- A -_ D - A - ] slower 
with increasing \U e «\ 

| i/off | > E a i single-electron transfer 

k[D- A- -►£>- A 2 -} 

and k [D 2- A -^ D - A -] 
slower with increasing \U e g\ 

\U e s\ > A, \U e s\ > Eal, 

T < | [/ e ff | electron-pair transfer 



TABLE I: Summary of the results. The effective Coulomb re- 
pulsion is denned by [/eff = U — V — 2au c . The corresponding 
reorganization energy is E a i = 2au> c the bias is e = 0. Start- 
ing with two electrons on the donor the system performs ei- 
ther a sequential single-electron transfer (D 2 ~ A — » D~ A~ — > 
DA 2 ~) or a pair transfer (D 2 ~ A — * DA 2 ~) depending on U e s ■ 

common bosonic bath. These interaction parameters can 
be calculated by ab initio methods for a specific system 
(see, for example ) 22 i 24 ). In our two-site model only the 
difference U = U — V enters the dynamics. In the pres- 
ence of a bosonic bath, the effective energy U is renor- 
malized to U Q s = U — 2aoj c . An effective attractive inter- 
action U c g < favors the localization of two electrons on 
the same site, a repulsive U c g > favors the distribution 
of electrons on different sites. 

The intricate correlated dynamics of two electrons de- 
pends on the activation energy. This is because, for 
the tunneling of an electron between two states, energy 
fluctuations are necessary for the reorganization of the 
donor-acceptor system and is influenced by an energy 
difference between the states. Therefore, the transfer 
characteristics in the unbiased case depends strongly on 
the effective on-site Coulomb repulsion U e g. Three rates 
have to be considered: the forward and backward rate 
between the double occupied states (D 2 ~A, DA 2 ~) and 
the two intermediate degenerate states (D~ A~) as well 



14 



as the direct rate between D 2 ~ A and DA 2 ~ . How these 
rates depend of U g is summarized in Table I. 

We have performed calculations for the probabilities 
P(t) of doubly and singly occupied donor and acceptor 
states using the time-dependent Numerical Renormal- 
ization Group metho d 33 ' . This information helps us 
to identify conditions under which the systems performs 
(a) concerted two-electron transfer, (b) uncorrelated, se- 
quential single-electron transfer or (c) fast concerted two 
electron followed by a single-electron transfer. With 
the time-dependent NRG method we can describe the 
crossover from damped coherent oscillations to incoher- 
ent relaxation as well as to localization (at T — > 0). The 
temperatures are chosen to be 0.1cj c > T > 3 • 10~ 8 u: c . 
For larger temperatures, when the bosonic bath can be 
treated classically, the Marcus rates are applicable. 

For U c s 3> A,E a i,T concerted electron transfer oc- 
curs in both methods: in the nuclear tunneling regime 
within the NRG as well as in the limit of a classical bath 
within the Marcus theory. As long as T < U e g, how- 
ever, thermal activation is absent and nuclear tunneling 
is the main process. Only a full quantum mechanical 
calculation yields the correct relaxation rates which are 
governed by quantum-fluctuation, dephasing and energy 
exchange with the environment. 

For small A/|{/ C ff| we found an effective pair hopping 
via virtual population of the low lying or high lying states 
D~A~. When the equilibrium probability for the states 
D~ A~ is finite, a slow single-electron accompanies the 
faster pair transfer. In contrast to the single-electron 
transfer with a frequency of the order A, the frequency 
of the pair transfer is of the order 4A 2 /|[/ C ff |. 

The concerted transfer becomes more uncorrelated and 
sequential at short times at high temperatures (T > 
U c tf), increasing coupling to the bosonic bath (E a i > 
U c s) or larger single-electron hopping (A > U e e)- The 
sequential transfer rate is non-monotonic with increasing 
U c tf- At first, the transition rate from D 2 ~ A to the delo- 
calized states D~A~ increases for small U e s > 0, reaches 
a maximum for U e s — E a \ before it decreases again. 
The rate for the consecutive process D~A~ — > DA 2 ~, 
however, decreases with increasing U. For a negative ef- 



fective Coulomb matrix element U c s, the transfer rate 
of the second process D~A~ — > DA 2 ~ is maximal for 
U e s = —E a i. In this parameter regime we expect that 
the second electron follows very shortly after the first 
electron was transferred. 

The transfer kinetics of more than two excess charges 
in, for example, biochemical reaction schemes or molecu- 
lar electronics applications is controlled by the molecule 
specific Coulomb interaction and its polar environment. 
Our study reveals the conditions for concerted two- 
electron transfer and sequential single-electron trans- 
fer. Concerted two-electron transfer is expected in com- 
pounds where the difference of the inter-site Coulomb 
repulsion and effective on-site repulsion are much larger 
than the single-electron hopping and larger than the 
temperature and reorganization energy. Furthermore, 
we have shown that non-monotonic characteristics of se- 
quential single-electron transfer strongly depends on the 
Coulomb interaction. A further study will include the in- 
fluence of a finite energy difference e between the donor 
and acceptor site. We will also report on the influence 
of Coulomb repulsion and many-particle effects on the 
long-range charge transfer using a longer Hubbard chain 
as bridge between donor and acceptor centers. 
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